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ABSTRACT 

Anomalies  in  surface  wave  patterns  observed  on 
SEASAT  Synthetic-Aperture  Radar  (SAR)  imagery  over 
Cobb  Seamount  are  examined.  These  anomalies  are  the 
surface  expression  of  currents  interacting  with  the 
seamount.  The  wavelengths  of  the  imaged  surface 
patterns  are  extracted  with  advanced  SAR  digital  image 
processing  techniques.  Nine  of  the  ten  detected 
wavelengths  match  the  expected  spectral  wavelength 
envelope  of  the  internal  waves  that  are  calculated 
using  oceanographic  data.  The  results  show  that  the 
wavelengths  of  the  surface  pattern  imaged  are  actually 
the  surface  manifestation  of  internal  waves  on  a 
subsurface  density  interface.  The  internal  waves  are 


created  by  current  flow  over  and  around  Cobb  Seamount 
disturbing  the  density  interface. 
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INTRODUCTION 


I.  h^k  ■  .Tfc  ^  W  M  UTh  W  « 


Subsurface  ocean  currents  flowing  over  underwater 
features  have  been  shown,  by  Lee  and  Beardsley  (1974) 
among  others,  to  alter  the  surface  wave  pattern.  The 
alteration  of  the  surface  wave  pattern  is  caused  by 
the  underwater  feature  disrupting  the  current  flow 
which  translates  by  hydrodynamic  processes  to  the 
surface  of  the  ocean  with  consequential  alteration  of 
the  surface  waves.  Satellite  observations  of  the  wave 
alterations  using  a  synthetic-aperture  radar  sensor 
have  been  shown  by  Kasischke  et.  al.(.1984)  to  have  a 
high  correlation  to  the  underlying  bathymetry 
(Fig.  1).  Similar  results  have  also  been  reported  by 
other  scientists  in  many  shallow  water  areas  with 
currents  (Fu  and  Holt,  1982;  McCandless  and  Mrazek, 
1982). 


The  alteration  of  surface  wave  patterns  in  the 
vicinity  of  deep  water  features  has  been  reported  by 
Robinson  (1985),  but  the  literature  does  not  report  a 
direct  correlation  between  imaged  patterns  and 
bathymetry  as  can  be  seen  in  the  imagery  of  Rockaway 
and  Ormonde  Seamounts  (Fig. 2  &  3).  Assuming  that  the 
same  hydrodynamic  dependence  on  subsurface  currents 
that  applies  to  shallow  water  imaging  applies  to  deep 
water  imaging,  a  research  project  was  undertaken  to 


Figure  2 

Rockaway  Seamount  (800  meter  depth) 
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examine  the  surface  patterns  associated  with  current 
flow  over  Cobb  Seamount  and  relate  the  surface 
patterns  to  specific  hydrodynamic  phenomena. 

HISTORY  AND  BACKGROUND 

Dr.  Gifford  Ewing,  an  oceanographer  well  known  for 
his  creative  thinking,  in  the  early  60's  postulated 
the  use  of  remote  platforms  in  the  detection  and 
measurement  of  oceanographic  conditions  and 
data.  Ewing  (1965)  noted  that  until  this  time  the 
thought  of 

"oceanography  from  a  satellite  sound(ed) 
incongruous  and,  to  a  generation  of 
scientists  accustomed  to  Nansen  bottles 
and  reversing  thermometers,  absurd". 

This  was  because  most  oceanographic  observations  until 

then  were  performed  through  the  laborious  process  of 

using  in-situ  instruments  and  attempting  to 

extrapolate  from  the  local  changes  the  larger  scale 

phenomena  creating  the  changes.  However,  with  the 

introduction  of  remote  sensors  the  possibility  of 

observing  synoptic  scale  phenomena  was  created. 

The  first  remote  sensors  were  visible  wavelength 
sensors  (such  as  cameras)  used  aboard  aircraft. 

Visible  wavelength  sensors  are  useful  in  the  detection 
of  current  fronts,  wave  heights,  and  wind  directions 
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(Duntley,  1965).  Visible  wavelengths  are  very 
susceptible  to  attenuation  in  the  atmosphere  and 
remote  sensing  is  possible  only  if  the  area  to  be 
imaged  is  naturally  illuminated  by  the  sun  or  moon  at 
the  time  of  imaging. 

Infrared  wavelength  sensors  (familiar  to  many  from 
weather  satellite  imagery),  are  one  of  the  most 
commonly  used  sensor  types.  Infrared  (IR)  wavelength 
sensors  (1  -  10  micrometer  wavelength)  are  used  to 
measure  relative  differences  in  thermal  activity  on 
the  Earth's  surface  and  in  the  atmosphere.  This 
allows  the  IR  sensors  to  detect  clouds,  surface 
currents,  warm  water  and  cold  water  interfaces,  and 
urbanized  areas.  Areas  of  colder  water  which  usually 
contain  higher  biological  activity  can  be  located  in 
this  manner.  IR  sensors  are  also  used  in  the 
development  of  daily  weather  forecasts  due  to  the 
ability  of  IR  sensors  to  image  cloud  formations  and 
fronts.  IR  sensors  are  not  dependant  on  direct 
illumination  from  the  sun,  as  visible  sensors 
are,  since  they  measure  the  amount  of  thermal 
radiation  emitted  by  the  surface  rather  than  reflected 
radiation.  However,  most  IR  wavelengths  are  highly 
attenuated  by  the  atmosphere,  particularly  by 
clouds.  This  attenuation  limits  the  capability  of 
these  sensors  in  imaging  the  oceans  in  areas  where 


storms,  fog,  or  cloud  cover  are  present. 

A  family  of  sensors  that  are  least  attenuated  by 
the  atmosphere  are  sensors  using  the  microwave 
wavelength  region  (wavelengths  of  0.1  to  100 
cm.).  The  longer  wavelengths  used  by  the  microwave 
sensors  are  only  scattered  by  the  very  large  rain 
drops  present  in  major  storm  clouds.  This  allows 
these  sensors  to  image  through  all  but  the  densest 
clouds.  The  first  space-borne  microwave  sensor  was 
the  radiometer  on  board  the  Mariner  2  that  was  used  to 
measure  surface  temperature,  atmospheric  water  vapor, 
ice  cover,  and  sea  temperature  on  Venus  (Carver, 
Elachi,  and  Ulaby,  1985).  One  type  of  centimeter 
wavelength  device  is  the  active  microwave 
sensor.  Those  sensors  do  not  require  outside 
illumination  because  the  sensors  send  out  energy  and 
measure  the  reflected  signal  energy. 

One  active  microwave  sensor  is  the  synthetic- 
aperture  radar  (SAR).  SAR  is  an  imaging  radar  that 
creates  a  very  high  resolution  (typically  10-20  m.) 
image  through  the  data  collection  and  processing 
technique  known  as  aperture  synthesis  (Appendix  A). 

Carl  Wiley  first  theorized  aperture  synthesis  as  an 
imaging  method  in  1951  while  working  for  Goodyear.  In 
1961  the  first  synthetic-  aperture  radar  sensor  was 
used  by  L.J.  Cutrona  to  image  ground  features 


(Harger,  1970).  A  space-borne  SAR  sensor  was  first 
used  to  image  the  oceans  in  1978  when  the  SEASAT 
satellite  was  launched.  Though  the  satellite  failed 
after  110  days,  100  million  square  kilome«.ers  of  SAR 
data  were  obtained  (Pravdo  et.al.,  1983). 

Included  in  the  SEASAT  data  are  the  imagery 
containing  surface  manifestations  of  the  underlying 
bathymetry.  The  information  drawn  from  the  imagery 
has  stimulated  researchers  to  attempt  to  model  the 
return  signal  created  by  current  flow  over  shallow 
water  bars  and  shoals  (Alpers  and  Henning,  1985).  The 
models  work  well  for  the  sensor  and  area  data  from 
which  they  were  derived,  but  none  has  proven 
successful  for  all  radar  wavelengths  and  conditions. 

SAR  imagery  in  areas  with  deep  water  topographic 
features  also  shows  surface  patterns  associated  with 
the  disruption  of  currents.  However,  since  the 
features  creating  the  disruption  are  deeper  than 
shoals  or  sand  bars  the  imagery  does  not  show  patterns 
conforming  to  the  underlying  bathymetry. 

A  commonly  imaged  surface  pattern  imaged  in  the 
region  of  deeper  water  features  is  an  internal  wave 
field  (Robinson,  1985).  For  example,  internal  waves 
are  often  created  by  current  flow  over  seamounts. 
Internal  waves  are  large  amplitude  waves  created  by 
the  disturbance  of  a  stable  density  difference 


interface.  Such  a  disturbance  could  be  caused  by 


current  flow  over  a  seamount.  Ewing  showed  that 
internal  waves  can  manifest  themselves  in  the  surface 
wave  pattern  (Fu  and  Holt,  1982).  These 
manifestations  of  internal  waves  in  the  surface  wave 
pattern  have  been  imaged  by  the  SEASAT  SAR  sensor 
(Robinson ,  1985). 


METHODS 


The  ideal  area  for  studying  the  surface 
manifestations  of  the  subsurface  current-seamount 
interaction  with  SAR  imagery  would  be  a  deep  flat 
ocean  plain  with  one  seamount  peak  reaching  near  the 
surface.  The  currents  should  be  steady  with  a  current 
meter  situated  on  the  seamount.  The  density  structure 
of  the  water  column  above  the  seamount  should  include 
a  strong  thermocline  (abrupt  change  in  the  temperature 
of  the  water  column  with  depth)  to  enhance  the 
potential  for  internal  waves.  The  surface  conditions 
should  include  a  stable  base  wave  pattern,  in  other 
words,  the  area  should  not  be  a  wave  generation  area 
since  the  surface  pattern  would  contain  too  many 
longer  wavelength  surface  waves  with  no  prevalent 
direction  of  propagation  for  analysis.  This  means 
that  the  ideal  wind  conditions  would  be  lower,  with 
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the  winds  only  strong  enough  to  develop  the  necessary 
capillary  waves.  Capillary  waves  are  short  wavelength 
waves  (less  than  2  cm.).  Capillary  waves  are  critical 
for  SAR  imaging  because  they  are  the  waves  that 
backscatter  the  transmitted  pulse  back  to  the  sensor 
to  create  the  image.  Of  course,  the  presence  of 
multiple  SAR  images  of  the  region  is  also  important. 
Since  not  all  SAR  data  shows  the  high  degree  of  detail 
that  would  be  preferable  for  analysis,  it  is  best  to 
find  an  area  that  has  as  many  passes  with  strong 
details  visible  as  possible. 

Cobb  Seamount  was  chosen  for  this  project  since  it 
most  closely  matches  the  ideal  conditions  set  forth. 
The  seamount  peak  protrudes  from  the  edge  of  the 
Cascadia  Abyssal  Plain  off  the  Oregon  coast  to  within 
30  meters  of  the  surface.  There  are  no  seamounts  in 
the  direct  path  of  the  North  Pacific  Drift  before  the 
current  interacts  with  Cobb  (Fig  4).  The  seamount  is 
located  in  the  North  Pacific  Drift  which  provides  a 
steady  1/3  Knot  of  current  in  the  month  of  July 
(Maughan,  1963).  The  tidal  currents  have  not  been 
directly  measured,  but  the  tidal  height  which  is  the 
actual  wave  created  by  the  tidal  forces  has  been 
measured  on  Cobb.  By  assuming  a  zero  time  delay 
between  the  tidal  height  and  the  tidal  current 
(reasonable  in  the  open  ocean),  an  estimate  of  the 
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cidal  contribution  to  the  current  can  be  found.  This 


is  estimated  to  be  up  to  one  knot  vectorially  added  to 
the  steady  N.  Pacific  Drift  current  (Apel,  1986). 

The  density  distribution  over  the  seamount  was  not 
directly  available  for  Cobb  at  the  time  of  SAR 
imaging,  but  temperature  readings  were  available  from 
Fleet  Numerical  Oceanography  Center  (FNOC),  Monterey, 
California.  This  data  shows  a  drop  of  7.5  degrees 
Celsius  between  21  and  49  meters  in  depth,  where  the 
peak  of  the  seamount  is  located.  The  data  was 
compared  to  seasonal  buoyancy  frequency  tables 
compiled  by  Bell,  Mays,  and  DeWitt  (1974).  The 
buoyancy  frequency  tables  (Fig.  5)  show  the 
Brunt-Vaisala  frequency  distribution  by  month  and  by 
region.  The  Brunt-Vaisala  frequency  is  the  natural 
harmonic  frequency  which  internal  waves  and  other 
density  instability  waves  resonate.  The  temperature 
data  showed  similar  changes  in  the  water  density  by 
depth  as  the  seasonal  data.  The  temperature  data 
available  only  contained  readings  at  large  steps  (over 
10  meters),  and  only  to  limited  depths  (less  than 
100  meters),  so  the  seasonal  data  were  used  in  the 
ca lculat  ions  . 

The  SAR  data  used  for  this  study  was  from  the  SAR 
sensor  on  board  the  SEASAT  satellite  launched  in 
1978.  Using  a  geo-location  program  from 


Winter 


Spring 


Summer 


Figure  5 

Seasonal  buoyancy  frequency  tables. 
(From  Bell,  Mays,  and  DeWitt,  1974) 
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Jet  Propulsion  Laboratory  (JPL),  it  was  determined 
that  only  three  recorded  SAR  passes  included  data 
taken  over  Cobb  Seamount.  The  three  data  sets  were 
included  in  satellite  revolutions  222,  466,  710.  The 
data  were  originally  transmitted  from  the  satellite  to 
the  Goldstone,  California  ground  station  where  the 
data  were  recorded. 

From  the  recorded  computer  tapes  of  this  data 
optical  and  digital  image  processing  can  be  performed 
(Fig.  6).  The  data  from  all  passes  is  available  on 
optical  film  for  the  SEASAT  SAR.  From  inspection  of 
the  three  revolutions  it  was  decided  that  revolution 
222  contained  too  little  detail  for  the  image  to  be  of 
use  in  the  analysis.  The  data  from  over  Cobb 
contained  in  revolutions  466  and  710  appeared  on  the 
optical  film  to  contain  sufficient  detail  for 
analysis.  The  data  from  466  and  710  were  processed 
digitally  and  prepared  for  analysis. 

Digital  processing  involves  taking  the  original 
transmitted  data  and  completing  the  entire  coherent 
addition  process  for  the  return  pulses  that  contain 
data  from  over  Cobb  Seamount.  The  output  image  is 
then  a  digital  image  with  each  resolution  element 
having  its  own  identity  and  intensity.  A  linear 
enhancement  was  also  performed  to  improve  contrast 
within  the  image.  The  digital  processing  was 
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Figure  6 

Optically  and  digitally  processed  imagery  of  SAR  data. 
(Imagery  from  Fu  and  Holt,  1982) 
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performed  on  the  SAR  digital  processing  system  at 
JPL,  with  the  enhancements  being  performed  at  Naval 
Research  Laboratory  /  Digital  Image  Processing 
Laboratory  (NRL/DIPL),  Washington,  D.C.. 

The  surface  conditions  present  when  the  Rev.  466 
data  were  taken  was  favorable  for  SAR  imaging.  The 
winds  at  the  time  of  imaging  were  light  according  to 
data  from  FNOC.  A  well  developed  wave  pattern  was 
already  present  in  the  area.  For  the  710  data, 
however,  the  winds  were  18  knots.  This  could  have 
caused  some  of  the  detail  in  the  surface  pattern  to 
have  been  obscured.  Because  of  this  obscuration,  it 
was  decided  to  drop  710  from  further  analysis  and 
focus  on  the  imagery  from  the  466  data  (Fig.  7) 

To  analyze  the  interaction  of  the  current  with  the 
seamount  and  the  surface  manifestation  of  the 
interaction,  accurate  bathymetric  data  was  obtained. 
Accurate  bathymetry  for  Cobb  Seamount  was  available 
from  the  Naval  Oceanographic  Command,  Bay  St.  Louis, 
Mississippi.  This  data  allowed  the  determination  of 
exact  location  of  the  seamount  and  aided  in  modeling 
the  interaction. 

A  summary  of  the  data  collected  for  this  project  is 
contained  in  Table  1. 
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DATA  SUMMARY 

DATA:  VALUE/FORM:  SOURCE  OR  DERIVED  FROM 

Bathymetry  Charts  Cobb,  Rockaway,  and  Ormonde 

Seamount  charts  obtained 
from  Naval  Oceanographic 
Office,  Bay  St.  Louis,  MS. 

Steady  1/3  Knot  Maughan,  1963. 

C  urrent 

Tidal  2/3  Knot  From  an  estimated  maximum 

Current  current  of  1  knot.  The 

tidal  height  is  two  hours 
from  max.  height 
(Schwinderski,  1981).  The 
estimate  was  developed  with 
the  aid  of  Dr.  Apel  (APL). 

Temperature  XBT  Data  From  Fleet  Numerical 

Distribution  Oceanographic  Center, 

Monterey,  CA. 

Brunt-  N  From  Bell,  Mays,  and 

Vaisala  DeWitt,  1974. 


466  SAR 
Data 


Original  data  from  Jet 
Propulsion  Laboratory, 
Pasadena,  CA. 

Data  enhanced  by  NRL/DIPL, 
Washington,  D.C.. 


Power  NRL/DIPL,  Washington,  D.C.. 

Spectrum 


Table  1 


Summary  of  collected  data. 


RESULTS 


From  examining  imagery  of  Rockaway  and  Ormonde 
Seamounts  (Fig.  2  &  3),  a  qualitative  model  was 
developed  for  comparison  to  SAR  imagery  taken  over 
Cobb  Seamount.  The  model  was  developed  to  aid  in  the 
understanding  of  the  dynamics  involved  in  the 
interaction  between  the  current  and  the  seamount,  and 
the  manifestation  of  the  interaction  in  the  surface 
wave  pattern.  The  model  assumes  a  homogeneous 
background  wave  pattern  onto  which  the  manifestation 
will  be  imprinted.  The  manifestation  is  expected  to 
have  two  major  constituents:  bow  waves  and  internal 
waves  (Fig.  8). 

Bow  waves  are  a  s em i pa r abo 1 i c  wave  pattern  created 
by  the  slowing  of  a  current  front  by  an  obscuring 
feature,  and  the  flow  of  the  current  around  the 
created  lower  velocity  area.  The  waves  form,  in 
the  case  of  a  point  stoppage  of  the  current  as  would 
be  created  by  a  seamount,  as  semiparabolic  waves 
pointing  into  the  current  and  cupping  around  the 
feature.  Since  Cobb  reaches  to  within  30  meters  of 
the  surface,  a  well  developed  bow  wave  pattern  would 
be  expected  to  form.  Such  waves  should  be  visible  in 
the  SAR  imagery  since  the  change  in  current  velocity 
and  the  associated  current  shear  will  affect  the 
surface  wave  pattern.  The  pattern,  however,  can  best 


SURFACE  PATTERN  MODEL  -  SIDE  VIEW 


SURFACE  PATTERN  MODEL  -  TOP  VIEW 


be  analyzed  through  qualitative  analysis  since  the 
bow  waves  do  not  have  a  constant  wavelength  or  common 
direction  of  propagation.  The  bow  waves  will  not 
propagate,  as  surface  waves  would,  since  they  are  only 
the  manifestation  of  a  current  shear.  If  the  waves 
did  have  a  constant  wavelength  and  common  direction  of 
propagation,  they  could  be  quantitatively  described  by 
using  an  analytical  methodology.  A  commonly  used 
methodology  is  to  develop  a  power  spectrum  of  the 
surface  wave  pattern.  A  power  spectrum  is  a  transform 
of  the  digital  information  into  a  spectrum  of  energy 
returned  versus  direction  of  propagation  and 
wavelength  of  the  surface  patterns  imaged. 

The  possibility  of  internal  waves  is  present  in  the 
area  of  Cobb.  Internal  waves  are  created  when  a 
strong  density  discontinuity  (such  as  a  strong 
thermocline)  exists  and  the  interface  of  the 
discontinuity  is  disturbed.  A  strong  thermocline  is 
present  in  the  area  of  Cobb  Seamount  and  the  current 
flowing  over  Cobb  would  disturb  the  interface. 

Robinson  (1985)  has  stated  that  internal  waves  have 
often  been  imaged  by  SAR  sensors  in  the  vicinity  of 
seamounts.  In  the  imagery  taken  of  Rockawav  Seamount, 
with  the  seamount  peak  located  800  meters  below  the 
ocean  surface,  internal  waves  are  in  evidence 
(Fig.  2).  Internal  waves  are  also  in  evidence  in  the 


imagery  of  Oromonde  Seamount  (25  meter  depth) 

(Fig.  3).  Although  Cobb  is  not  as  shallow  as  Ormonde 
or  in  the  very  strong  flow  of  the  Gulf  Stream  like 
Rockaway,  it  was  predicted  that  internal  waves  would 
be  imaged  over  Cobb  also. 

Current  interaction  with  the  normal  development  of 
internal  waves,  however,  affects  the  internal  wave 
wavelengths  that  will  manifest  themselves  in  the 
surface  wave  pattern.  Internal  waves  created  with 
minimal  current  flow  or  by  simple  displacement  form 
with  wavelengths  determined  by  the  natural  buoyancy 
frequency  ( B r un t-Va i sa la )  determined  by  the  density 
changes  present.  If  a  strong  current  is  present,  the 
internal  waves  are  forced  to  oscillate  at  higher 
frequencies.  This  forcing  causes  the  internal  waves 
to  form  in  an  envelope  of  frequencies  in  the  lee  of 
the  seamount.  The  envelope  has  a  lower  wavelength 
cutoff  determined  by: 

Wavelength  = _ C _  E q  .  1 

N 

where  C  is  the  magnitude  of  the  current 

N  is  the  effective  B r un t - Va i s a  1  a  frequency 

(Apel, 1986) . 

The  other  end  of  the  wavelength  envelope  is  the 
determined  from  the  B r un t - Va  i  s a  1  a  frequency 
alone.  The  internal  waves  in  this  envelope  will  form 
in  the  lee  of  the  seamount,  then  when  the  current 
abates  the  waves  will  travel  over  the  seamount  and 


propagate  out  into  the  ocean  in  the  direction 
opposite  the  original  creating  current  (Apel, 

1986).  Using  seasonal  data  for  N  in  Eq.l  and  a 
current  estimate  of  1/3  knot  steady  current  and  2/3 
knot  tidal  current  adding,  the  lower  end  cutoff 
wavelength  is  0.5  km.  The  imagery  is  therefore 
expected  to  contain  manifestations  of  internal  waves 
of  variable  wavelengths  with  all  the  wavelengths 
greater  than  one  half  a  kilometer. 

The  combination  of  a  bow  wave  pattern  and  the 
surface  manifestation  of  the  internal  wave  generation 
area  is  taken  as  the  qualitative  model  of  the  expected 
imagery.  A  possible  version  of  the  model  is  shown  by 
figure  8  . 

The  revolution  466  digital  imagery  does  show  some 
similarities  to  the  proposed  model.  There  is  a  subtle 
signature  from  two  bow  waves  (Fig.  9),  although  not 
nearly  as  clearly  defined  a  signature  as  was 
expected.  There  is  also  an  area  that  shows  the 
presence  of  longer  wavelength  manifestations  inside 
the  bow  waves.  The  pattern  is  confused  as  would  be 
expected  from  the  combination  of  multiple  wavelength 
waves  in  one  region  (Fig.  10). 

Performing  a  power  spectrum  analysis  on  the 
multiple  wavelength  waves  is  a  proven  method  of 
quantifying  the  patterns  imaged.  The  power  spectrum 


Bow  wave  manifestations  on  466  enhanced  imagery. 
(Dark  line  is  drawn  between  the  two  manifestations 
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most  commonly  used  is  a  representation  of  the 
returned  power  versus  the  wavelength  and  direction  of 
propagation  of  the  waves  returning  the  power.  The 
power  spectrum  used  for  this  analysis  were  generated 
through  the  use  of  a  Fast  Fourier  Transform  contained 
in  the  software  in  the  mainframe  computers  at  Naval 
Research  Laboratory  (NRL/DIPL). 

A  power  spectrum  was  developed  of  two  areas  in  the 
466  image  (Fig.  11):  a  reference  area  that  contains  no 
manifestations  of  the  current-seamount  interaction, 
and  a  test  area  that  is  centered  in  the  manifestations 
of  the  interaction  (Fig.  12  &  13).  From  the  power 
spectrum  it  was  determined  that  the  predominate 
direction  of  wave  propagation  was  to  or  from  the 
North-Northeast.  Using  this  orientation,  a 
directional  slice  was  taken  through  the  power  spectrum 
(Fig.  14)  to  determine  a  two  dimensional  power 
spectrum  (power  vs.  wavelength)  (Fig.  15  &  16).  These 
represent  both  the  positive  and  negative  directions  of 
propagation  possible  in  the  image  (180  degrees 
opposite).  Since  there  should  be  perfect  symmetry 
about  the  spectrum  centerpoint  for  waves,  one  half  of 
the  spectrum  showing  only  the  positive  direction 
wavelengths  can  be  plotted,  and  expanded 
(Fig.  17  6  18).  These  can  then  be  combined  for  easy 
comparison  (Fig.  19). 
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Figure  15 


Full  power  spectrum  of  reference  area. 
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POWER  SPECTRUM-COBB  SEAMOUNT  AREA 
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POWER  SPECTRUM-COBB  SEAMOUNT  AREA 
SEASAT  SAR  REV  466  IMAGE 
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Figure  18 

Half  power  spectrum  of  test  area. 
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EXPANDED  POWER  SPECTRUM  -  COBB  SEAMOUNT  8  REFERENCE  AREA 
SEAS AT  SAR  REV  466  IMAGE 
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Figure  19 

Comparison  of  half  power  spectrums  of  reference 
and  test  areas. 
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From  the  combined  half  power  spectrums  it  can  be 
shown  that  there  is  a  significantly  higher 
contribution  to  the  power  spectrum  from  longer 
wavelength  waves  in  the  target  area  as  opposed  to  the 
reference  area  (Fig.  19).  In  particular,  there  are 
significant  (relative  gain  of  9  dB  or  more) 
differences  in  the  power  associated  with  peaks  at  ten 
different  wavelengths  (Table  2).  The  standard  cut  off 
value  for  determining  significance  is  a  6  dB  gain  over 
system  noise  (Mango,  1986).  The  system  noise  can  not 
be  determined  from  one  image,  but  the  noise  level 
should  be  the  same  for  both  areas.  The  9  dB  relative 
gain  level  was  selected  to  ensure  that  any  changes  in 
system  noise  between  the  two  image  areas  would  not 
effect  the  significance  of  peaks  in  the  spectrum. 

There  are  also  peaks,  with  less  than  9  dB  gain,  at 
other  long  wavelengths,  but  these  could  be  attributed 
to  noise  in  original  receipt  of  data,  downlink 
transmission,  and  processing.  There  are  no  long 
wavelength  significant  peaks  in  the  reference  area 
spec  t  rum . 

All  of  the  wavelengths,  except  the  0.48  km. 
wavelength,  are  enclosed  in  the  predicted  internal 
wave  wavelength  envelope.  The  0.48  km.  wavelength 
present  outside  the  predicted  envelope  is  most  likely 
due  to  an  error  in  the  current  estimate,  placing  the 


SIGNIFICANT 

WAVELENGTHS 

Wavelength 

(km.) 

Power 

Target  Area 

Power 

Reference  Area 

Relative 

Gain 

7.26 

12.52 

9.16 

14.6 

3.61 

11.59 

8.16 

14.9 

2.51 

10.70 

8.40 

10.0 

1.91 

9.70 

6.50 

13.9 

1.58 

'9.35 

6.78 

11.2 

0.89 

8.70 

6.40 

10.0 

0.85 

8.81 

6.00 

12.2 

0.60 

7.40 

5.20 

9.6 

0.55 

7.10 

4.90 

9.1 

0.48 

7.15 

4.80 

10.2 

Table  2 


Significant  wavelengths 


cut  off  wavelength  too  long  due  to  an  over¬ 


estimation.  The  wavelengths  are  much  too  long 
for  surface  wind  waves,  which  have  wavelengths  on  the 
order  of  tens  of  meters  or  less.  The  peaks  detected 
in  the  power  spectrum  must  be  from  the  surface 
manifestation  of  internal  waves  created  by  current 
flow  over  Cobb  Seamount. 

CONCLUSIONS 


Changes  in  the  surface  wave  pattern  are  imaged  by 
the  SEASAT  SAR  sensor  in  the  area  of  Cobb  Seamount. 

The  surface  changes  can  be  attributed  to  surface 
manifestations  of  hydrodynamic  phenomena  associated 
with  current  flow  over  seamounts.  The  particular 
phenomena  imaged  are  bow  waves  and  internal  waves. 
Subtle  manifestations  of  the  bow  waves  are  visible  on 
the  enhanced  image  indicating  that  they  are  present, 
but  not  highly  developed  in  the  Cobb  Seamount  area. 

The  internal  waves  were  predicted  to  be  contained  in 
an  envelope  of  wavelengths  greater  than  0.5  km.  Large 
wavelength  surface  patterns  are  visible  in  the  imagery 
and  are  quantitatively  measured  using  a  power 
spectrum.  Nine  of  the  ten  significant  peaks  found  in 
the  test  area  power  spectrum  are  contained  within  the 
calculated  envelope.  The  presence  of  a  strong 


thermocline  and  current  in  the  region  of  Cobb 
Seamount  support  the  conclusion  that  the  peaks  in  the 
power  spectrum  are  the  surface  manifestations  of 
internal  waves  created  by  current  flow  over  Cobb 
Seamount . 


APPENDIX  A 


S y n t he t ic -a p e r t ur e  radar  (SAR)  is  a  high  resolution 
imaging  radar  system  which  includes  a  sensor,  data 
recording  system,  and  a  ground  image  processing 
computer  system.  The  spacial  resolutions  achieved  by 
SAR  systems  is  on  the  order  of  tens  of  meters  as 
opposed  to  the  kilometer  resolutions  of  standard  radar 
sensors.  The  resolution,  or  ability  to  distinguish 
closely  spaced  objects  in  an  image,  is  measured  in  two 
different  directions:  across  the  satellite  ground 
track  (range)  and  along  the  satellite  track  (azimuth). 

The  across  track  or  range  resolution,  in  the  case 
of  the  SEASAT  satellite  SAR,  is  improved  by  the 
particular  geometry  used  by  the  satellite  system 
(Fig.  Al).  Standard  radars  image  by  transmitting  an 
electro-magnetic  pulse  directly,  out  and  measuring  the 
reflected  energy.  The  SEASAT  SAR  instead  uses  a 
side-looking  geometry.  By  looking  to  the  side  of  the 
satellite  ground  track,  the  system  develops  a  time 
difference  between  returns  from  the  ground  imaging 
swath  regions  closer  to  the  satellite  track  (Fig.A2). 
The  resolution  is  determined  by: 

IFOVr  =  \_  _Ljd_*_c_  Eq.Al 

2  cos'©1 

where  IFOVr  is  the  instantaneous  field  of  view 
(resolution)  in  the  range  direction 
Lp  is  the  pulse  length  in  seconds 
c  is  the  speed  of  light 
■©■  is  the  grazing  angle. 


Satellite 
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The  SEASAT  SAR  system  has  an  across  track  resolution 
of  approximately  21  meters. 

Synthetic-aperture  radar  differs  from  traditional 
radars  principally  in  the  along  track  resolution.  The 
equation  governing  radar  resolution  for  the  along 
track  (azimuth)  direction  is: 

IFOVaz  -  J.  X*  R  Eq.A2 

2  Daz 

where  IFOVaz  is  the  instantaneous  field  of  view 
in  the  azimuth  direction 
A  is  the  wavelength  transmitted 
R  is  the  slant  range  from  the  sensor  to 

the  ground  imaging  swath  (Fig.  A2) 
Daz  is  the  antenna  diameter  in  the 
azimuth  direction. 

For  the  SEASAT  SAR  sensor  this  would  result  in  a 
resolution  of  9  kilometers.  The  resolution  is 
improved  to  the  tens  of  meters  desired  by  creating  an 
antenna  much  larger  then  the  11m.  by  2m.  physical 
antenna  carried  by  the  SEASAT  SAR  sensor  (Fig.  A3). 
This  antenna  is  an  artificially  created  antenna  or 
synthetic-aperture  . 

S yn t he t i c- ap e r t u r e s  are  created  by  coherently 
adding  multiple  return  pulses  from  an  imaged  area  as 
the  satellite  travels  through  space.  When  the 
multiple  pulses  are  coherently  added  together  the 
processing  system  creates  a  single  image.  The  image 
has  a  resolution  equal  to  the  resolution  that  would  be 
possible  using  an  antenna  the  size  of  the  synthetic- 
aperture.  The  diameter  of  the  s y n t h e t  i  c -ap e r t u r e  is 
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equal  to  the  distance  the  satellite  traveled  while 
transmitting  and  receiving  the  multiple  pulses.  The 
developed  synthetic-aperture  length  (Lsa)  is: 

L  s  a  =  ^  *  R  E  q  .  A  3 

Da  z 

where  Lsa  is  the  s y n t he t  ic -a p e r t u r e  length 
developed 

A.  is  the  wavelength  transmitted 

R  is  the  slant  range 

Daz  is  the  diameter  of  the  physical 
antenna. 

The  multiple  pulse  addition  that  creates  the 
synthetic-aperture  is  known  as  coherent  addition. 
Coherent  addition  is  performed  by  the  satellite 
retaining  the  Doppler  history  of  a  ground  target's 
multiple  returns  as  the  ground  target  travels  through 
the  ground  imaging  swath.  This  Doppler  history  is 
then  used  to  separate  the  signal  return  from  a  small 
area  (on  the  order  of  20  meters  across). 

The  Doppler  history  is  a  recording  of  the  multiple 
returns  from  a  single  area  as  both  frequencies  and 
intensities.  As  the  satellite  travels  through  space, 
a  relative  velocity  is  created  between  the  satellite 
and  the  earth's  surface  area  the  satellite  is 
imaging.  This  relative  velocity  creates  a  shift  in 
the  frequency  returned  by  the  area,  to  the 
satellite.  This  frequency  shift  is  known  as  a  Doppler 
shift.  If  the  area  to  be  imaged  is  ahead  of  the 
satellite  the  area  will  have  a  positive  shift  in 


frequency  as  it  enters  the  ground  imaging  swath.  As 
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the  satellite  continues  through  the  imaging  swath  the 
Doppler  shift  will  decrease  until  the  area  is 
perpendicular  to  the  satellite  track.  When  the  area 
is  perpendicular,  the  Doppler  shift  will  be  zero,  and 
the  return  signal  will  have  the  same  frequency  as  the 
original  transmitted  pulse.  The  Doppler  shift  will 
continue  decreasing  as  the  area  goes  through  the 
imaging  swath,  with  the  shift  being  an  increasingly 
negative  value  (Fig.  A4 ) . 

The  SAR  system  not  only  records  the  frequency 
Doppler  history  of  each  ground  area,  but  also  the 
return  signal  intensity  associated  with  each  return 
(Fig.  A5).  The  coherent  addition  process  takes  the 
return  signals,  including  both  frequency  and 
intensity,  and  adds  the  returns  to  determine  the 
intensity  of  return  when  the  Doppler  shift  equaled 
zero.  This  intensity  is  used  as  the  intensity  for  the 
pixel  area  in  the  image  produced. 

The  pixel  size  is  determined  by  the  resolution  of 


the 

system. 

Taking  the 

syn t he t i c -a pe r t u r e  length 

(Eq  . 

A3 )  as 

the 

antenna 

d i ame  ter  and 

substituting  that 

into 

Eq  .  A2 

the 

resulting  resolution 

is  found  : 

IFOVaz 

=  1 

A  *  R 

=  Da  z 

Eq  .  A4 

2 

L  s  a 

2 

For 

SEASAT 

this 

allows 

a  calculated 

resolution  of 

6  meters.  Usually  the  processing  system  does  not  take 
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A  point  on  the  ground  will  first  be  imaged  at  point  a  ^ 

The  ground  point  will  then  be  imaged  at  points  a  a  ^,et-c 
as  the  satellite  continues  on  its  velocity  track 


Figure  A4 

Doppler  shift  of  return  signal  (single  return  pulse). 


DOPPLER  HISTORY  OF  ONE  POINT  -  MULTIPLE  LOOK  (FREQUENCY) 


Figure  A5 

Doppler  history  of  the  multiple  returns  from  one  ground 
point  (frequency  and  intensity). 


full  advantage  of  this  achievable  resolution.  This 
is  done  to  limit  the  image  speckle  effect  intensified 
by  full  coherent  addition.  Instead  only  a  part  of  the 
full  synthetic-aperture  is  coherently  processed. 
Coherent  sub -ap e r t ur e s  are  then  added  together 
incoherently.  This  is  called  taking  looks  or 
multi-look  processing.  The  result  is  that  spacial 
resolution  is  compromised  to  reduce  speckle.  The 
images  in  this  report  were  four  look  processed  for  an 
along  track  spacial  resolution  of  about  24  meters. 
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